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Abstract —This paper investigates the virtualization and op¬ 
timization of a multi-cell WLAN. We consider the station 
(STA)-access point (AP) association and airtime control for 
virtualized 802.11 networks to provide service customization 
and fairness across multiple internet service providers (ISPs) 
sharing the common physical infrastructure and network ca¬ 
pacity. More specifically, an optimization problem is formulated 
on the STAs’ transmission prohahilities to maximize the overall 
network throughput, while providing airtime usage guarantees 
for the ISPs. Subsequently, an algorithm to reach the optimal 
solution is developed by applying monomial approximation and 
geometric programming iteratively. Based on the proposed three- 
dimensional Markov-chain model of the enhanced distributed 
channel access (EDCA) protocol, the detailed implementation 
of the optimal transmission probability is also discussed. The 
accuracy of the proposed Markov-chain model and the perfor¬ 
mance of the developed association and airtime control scheme 
are evaluated through numerical results. 

Index Terms —Wireless local area networks (WLANs), IEEE 
802.11e standard, STA-AP association, airtime control, through¬ 
put, fairness, Markov chain, geometric programming. 

I. Introduction 
A. Background and Motivation 

Wireless virtualization has recently emerged as an archi¬ 
tectural choice for the wireless networks, in which different 
service providers can share physical infrastructure and wireless 
resources. Allowing resources to be shared, vitualization can 
facilitate a flexible and lower-complexity solution to support 
customized services with finer control over quality-of-service 
(QoS) features. To enable service customization, the key issue 
would be to provide isolation among virtual networks running 
by different service providers. Such isolation can be achieved 
through contention-free resource allocation techniques based 
on TDMA and/or EDMA, by isolating resources across ser¬ 
vice providers. But, in systems using contention-based access 
protocols, providing isolation among virtual networks is more 
challenging. 

Eor instance, in virtualized 802.11 WLANs, transmissions 
of different virtual WLANs (V-WLANs) are closely coupled, 
although administrative virtualization (i.e., one physical AP 
advertises multiple service set identifiers (SSIDs)) can already 
differentiate groups of flows. With a carrier sense multiple 
access (CSMA)-based MAC, unavoidable collisions act to 
couple the transmissions of different V-WLANs. Moreover, 
since the network capacity is shared yet constrained, the 


increase of traffic in one V-WLAN may reduce the available 
network capacity to another ||T|. Thus, an efficient resource 
allocation among V-WLANs is essential to manage the MAC- 
layer couplings. 

In a 802.11 WLAN with densely deployed APs, before a 
STA can access the network, it needs to make a decision 
about which AP to associate with. In virtualized 802.11 
networks, such STA-AP association control could create an 
opportunity to provide fairness guarantees among different 
ISPs. In addition to STA-AP association, the airtime among 
STAs associated to the same AP can optimally be adjusted. 
Controlling airtime usage of the STAs provides the opportunity 
to optimize the WLAN performance (e.g., improving WLAN 
throughput by exploiting multi-user diversity) as well as 
another degree of freedom to guarantee fairness among the 
ISPs. 

The main focus of this work is on the virtualization and 
optimization of a multi-cell WLAN. More specifically, STA- 
AP association and airtime control are jointly explored to 
provide fairness and throughput guarantees for different V- 
WLANs. Taking into account STA transmission rates and ISP 
airtime reservations, an optimization problem is formulated to 
adjust the transmission probability of each STA at each AP. 
The objective is to maximize the overall network throughput, 
while keeping a total airtime guarantee for each ISP. To 
solve the formulated problem which is non-convex and thus 
computationally intractable, an iterative algorithm is developed 
through successive geometric programming. This algorithm 
can achieve an optimal solution with an affordable complexity. 

Eurthermore, to implement the optimal transmission proba¬ 
bilities, we model the EDCA protocol with a three-dimensional 
Markov chain and establish the relationship between the 
transmission probability of a STA and the EDCA param¬ 
eters. Based on this relationship, we developed a control 
algorithm to approach the optimal transmission probability by 
jointly manipulating the EDCA parameters such as contention 
window size (CW) and arbitration inter-frame space (AIES). 
Einally, through numerical results, we verify the performance 
enhancement provided by the developed STA-AP association 
and airtime control approach in terms of throughput and 
fairness. 


B. Related Works 

1) STA-AP Association and Airtime Control: In most cur¬ 
rent vendor implementations, 802.11 STAs choose the AP with 
the highest received signal-to-noise ratio (SNR) to connect 
with. Since the STA density is often uneven in the network 
0, 0, the Max-SNR approach can lead to an unbalanced 
distribution of STAs among APs, causing unfairness. In order 
to balance the load of APs, several STA-AP association 
algorithms have been presented in the literature, mostly by 
maximizing the minimum throughput of all STAs 0-0- Nev¬ 
ertheless, in a basic service set (BSS) including an AP and its 
associated multi-rate STAs, it is shown that the throughput is 
limited by the STA with the lowest data rate. This phenomenon 
is also to as the performance anomaly problem Thus, 
comparing with the Max-SNR approach, these load-balancing 
approaches improve the max-min fairness among STAs at the 
cost of decreasing the aggregate throughput. 

To address the performance anomaly and balance the trade¬ 
off between aggregate throughput and fairness, proportional 
fair throughput allocation has widely been considered in multi¬ 
rate 802.11 WLANs 0-|Tg. In g, proportional fairness is 
studied in a single BSS. It is shown that propotional fairness 
leads to an airtime-fairness, where equal airtime usage is 
provided to all STAs. Moreover, in a multi-AP WLAN, O). 
GD study STA-AP association problem with an objective to 
maximize the proportional fairness. More precisely, associa¬ 
tion control is implemented in a form of airtime allocation, 
where the transmission time of STAs at different APs are 
jointly optimized im^ im¬ 
precise controlling of the STA airtime usage in a 802.11 
contention-based WLAN is very difficult due to the distributed 
and random nature of their CSMA-based MAC protocol. Nev¬ 
ertheless, the airtime control algorithms have been proposed 
in the literature by manipulating the MAC parameters im- 
0- For instance, 0 discusses how to dynamically modify 
the CW based on the transmission data rate and estimated 
packet error rate (PER) in a WLAN aiming to guarantee 
airtime fairness. Similarly, in p^ , CW is adjusted according 
to the traffic load of each STA to improve throughput and 
fairness in the WLAN. Furthermore, 0 proposes a machine¬ 
learning approach to adjust AIFS and CW aiming to achieve 
airtime fairness and reduce frame delay. Such heuristic works 
provide pretty good understanding of the influence of specific 
MAC parameters such as CW and AIFS on the transmission 
probability and thus airtime usage of the STAs. However, the 
discussions are limited to adjusting the parameters separately 
and thus can only And a suboptimal solution in a subset of the 
feasibility region. 

In a virtualized WLAN serving multiple ISPs, STA-AP 
association and airtime control become more challenging. The 
reason is that fairness guarantees and service customization 
are required for ISPs, while there are unavoidable couplings 
among the STA transmissions of different ISPs in the network. 
There are a few works addressing only airtime control in the 
literature. Considering a virtualized single-cell WLAN, in 


a heuristic airtime control algorithm is proposed to achieve the 
target airtime usage for each ISP by controlling the minimum 
CW of each STA. Similarly, p^ addresses optimizing CW 
using control theory. But, as the discussion is limited to 
controlling minimum CW, the optimality of the result might 
be sacrificed. In a multi-cell WLAN, Q presents an analysis 
on the feasibility region of ISP airtimes and characterizes 
the ISP airtimes at the rate region boundary. Furthermore, a 
distributed algorithm is developed to allocate airtime slices 
among ISPs and flow rates within each slice in a max-min 
fair manner. Since max-min fairness is used as an objective 
for rate allocation among the flows in each ISP, the optimality 
of the achieved total throughput may not be guaranteed in a 
multi-rate WLAN. In addition, the association control is not 
discussed in ||T). 

2 ) IEEE EDCA Modeling and Optimization : To implement 
the optimal STA-AP association and airtime allocation, it is 
desirable to control the transmission probability of each STA. 
However, in a CSMA/CA-based WLAN, the only controllable 
parameters are the MAC layer parameters. Thus, it is essential 
to mathematically model the effects of such parameters on the 
transmission probabilities of STAs. 

The EDCA modeling has received considerable attention in 
the literature. There are several studies, such as 0, 0, 
that only focus on numerically solving the stationary state 
(transmission probability of each STA) given the network 
configuration. Another group of works, e.g., |[22)-|[29), provide 
models for the EDCA protocol with the help of a Markov 
chain, inspired by the classical work of Bianchi pO) . Using 
such Markov chain models, it is possible to optimize the 
network by tuning the EDCA parameters. However, due to 
the complexity of the proposed EDCA models, it is hard 
to establish an explicit relationship between the transmission 
probabilities and the controllable EDCA parameters. Thus, 
these works mostly provide a numerical (e.g., 0-0) or 
approximate analytical solution in terms of MAC parameters 
(e.g., |26|-||29l). Moreover, none of these works has jointly 


controlled the EDCA parameters. 

Our proposed Markov chain model is based on pT) , with 
which an explicit relationship between the transmission prob¬ 
ability and EDCA MAC parameters can be established. One 
limitation of the EDCA model proposed in pT) is the accuracy 
of AIFS differentiation. Unlike in DCF, the time duration each 
STA has to wait before its backoff process (i.e. AIFS) is 
different in EDCA. This can lead to the result that the number 
of contending STAs is not time-homogeneous pT| . Therefore, 
the collision probability that each STA faces is also not time- 
homogeneous. But, as pointed out in p7| , time-homogeneity 
assumption greatly simplifies the modeling complexity, and 
thus, the steady state performance can be characterized as 
explicit functions of backoff parameters. Furthermore, the 
accuracy of such model can be effectively improved by setting 
a relatively large initial CW. 

C. Structure 

The rest of this paper is organized as follows. Section |I^ 
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Fig. 1: Layered system model 



presents an overview of the system configuration and model¬ 
ing. In Section [I^ we first analyze the feasibility region of the 
transmission probabilities based on the proposed Markov chain 
model for IEEE EDCA. Then, we formulate the transmission 
probability optimization problem, which maximizes the system 
throughput and guarantees the fairness among the ISPs. In 
Section [TVl the implementation details of MAC parameter con¬ 
trol are discussed in order to achieve the optimal transmission 
probability. Illustrative results are provided in Section to 
evaluate the performance of the developed algorithms. Section 
|V^ provides some concluding remarks. 

II. System Configuration and Modeling 

We consider an IEEE 802.11-based WLAN that consists 
of a large number of APs. APs operate on non-overlapping 
frequency channels. Let A be the set of APs and Na = |,A| 
be the total number of APs. Each AP has a limited coverage 
area and all STAs are randomly distributed in the field. The 
network carries traffic belonging to a number of different ISPs 
(also referred to as V-WLANs). Let /C be the set of ISPs using 
the network. Eurthermore, let Sk be the set of STAs of ISP 
k G K, and Nk = |Nfc| be the number of STAs belong to ISP k. 
Eurthermore, let S be the set of all STAs and Ng = 
be the total number of STAs in the network. The network 
is administratively virtualized, i.e., each AP will broadcast 
multiple different SSIDs, one for each ISP. Pigure[^ illustrates 
an example of the network architecture with four physical APs 
and two ISPs. 
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Fig. 3: IEEE 802.lie EDCA channel access procedure 


A. Enhanced Distributed Channel Access (EDCA) 


To access the channel, the STAs are assumed to follow the 
EDCA MAC protocol of 802.1 le standard. EDCA protocol 
is designed to enhance the basic MAC mechanism, i.e., 
distributed coordination function (DCE), aiming to support 
service differentiation. Similar to DCE, EDCA is also a 
contention-based access scheme, based on CSMA/CA using 
binary exponential backoff rules to manage retransmission of 
collided packets. Here, we briefly review the EDCA operation. 


as standardized by 802.11 protocol |30|, 

EDCA requires a STA to monitor the channel before any 
transmission. If the channel is sensed idle for a time interval 
equal to an arbitration inter-frame space (AIES), the STA 
transmits. Otherwise, if the STA senses a transmission either 
immediately or during the AIES, it continues monitoring the 
channel. When the channel is measured idle for an AIES, 
the STA backoffs for a random period of time. EDCA uses 
a discrete-time backoff mechanism, i.e., the time following 
an AIES is slotted. The backoff time is selected according 
to a uniform distribution in the interval [0, Wi] where Wi 
represents the contention window of STA i S S |30|, 


STA i starts a packet transmission at backoff stage 0, 
where Wi is set equal to the minimum contention window 
size Wmin.i- Then, after each unsuccessful transmission, STA 
i moves to the next backoff stage and doubles Wi. When 
the maximum backoff stage is reached, Wi is no longer 
increased and stays at T^^Wjnia.i- If the STA experiences a 
collision at the backoff stage nii, it will retry transmission for 
at most hi times, where hi is the retransmission limit at the 
maximum backoff stage. If the packet is still not successfully 
transmitted after hi retransmissions, it will be discarded 1^ , 
il, (33). 


The backoff time counter is decremented and a STA trans¬ 
mits when the backoff counter reaches zero. If the data frame 
is successfully received, the AP waits for a period of time 
called short inter-frame space (SIPS) and then sends an ac¬ 
knowledgment (ACK). Once the channel is sensed busy while 
counting down, the STA will freeze its backoff counter and 
continue decrementing when channel finds idle again. Pigure 
illustrates an example of the channel-access procedure of 
two STAs using EDCA (3^, (^, (S^. 

EDCA allows to have variable waiting times and contention 
parameters (i.e., AIES and Wmin) to enable service differenti¬ 
ation. Using EDCA, for each transmission opportunity, a STA 
can transmit multiple back-to-back packets for a fixed period 



















of time. Let Ttxop be the duration of a data frame. Then, the 
duration of a successful transmission becomes 
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T, = Ttxop + SIFS + 7 + ACK + 7 + AIFS (1) 

where 7 denotes the transmission delay. Similarly, the duration 
of a collision can be calculated as 

T; = rTxop + 7 + AIFS. (2) 

Note that when the colliding STAs use different AIFS values, 
the AIFS value in should take the value of the largest 
AIFS. But, since the hxed transmission duration Ttxop (in the 
order of ms) is much larger than AIFS in the order of /rs, the 
difference between the AIFS can be ignored here. 


B. Association Control via Transmission Probabilities 


In a WLAN with APs densely deployed, STAs need to 
determine which APs to connect with. We aim to generalize 
the association control problem by adjusting the transmission 
probability of each STA at any AP, rather than selecting one 
AP to associate with. Thus, we define r“ (0 <rf< 1) as the 
probability that STA i attempts to transmit at AP a in a general 
time-slot. Consequently, the probability that a time-slot is idle 
in the BSS including AP a is 

In a given BSS, transmitted packets will be received success¬ 
fully, if exactly one STA transmits on the channel. Thus, the 
probability of a successful transmission initiated by STA i 
becomes 


■ = t‘ 
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(4) 


Since ACK and SIFS are relatively small (in the order of /rs) 
compared with Tjxop (in the order of ms), we approximate 
Tg and to be of the same size and denote them by 
T. Consequently, the expected length of a general time-slot 
becomes 


E{Tg} = SPS,,, + (1 - (5) 

where 5 is the duration of an idle time-slot. Furthermore, the 
expected information (in bits) transmitted by STA i to AP a 
in a general time-slot can be derived as 

E{4} = PLc,^rtTTXOP (6) 

where rf represents the transmission data rate of the link 
between STA i and AP a. As defined in H), based on the 
0 and (|^, the throughput of STA i at AP a becomes 

^ E{/g} ^ -Psucc,zL°^TXOP 

* E{rj P.^j + ii-P.^jT- 

Let define a new variable x“ = (x^ > 0), which 

represents the expected number of consecutive transmission 
attempts by STA i at AP a as g, 0. Consequently, P-^i^ 
and j will be transformed into 
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nPe5(l + <') 

Subsequently, from 0 and 0, T“ can be represented in terms 
of x“ as 
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where t = and t' = 

In addition to the throughput of each STA, the fraction of 
time that each STA spends for transmission could be consid¬ 
ered as another performance metric, specifically in order to 
measure and preserve fairness among different STAs or ISPs. 
The total access airtime for STA i-including both successful 
transmissions and collisions-becomes 
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where , = rf 1 - Tip 


is the probability 
that STA i suffers from a collision in a general time-slot as 
defined in Q. Consequently, 
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( 12 ) 


In this work, we aim to maximize the overall network 
throughput, while guaranteeing a minimum requirement on the 
aggregate airtime of each ISP. To this end, the transmission 
probability of STAs (t“) need to be adaptively optimized 
by maximizing the aggregate throughput of all STAs at all 
APs (i-e-,X]iG5 Furthermore, for each ISP (e.g., ISP 

k), a constraint needs to be set in order to keep the total 
airtime of all STAs belonging to ISP k larger than a minimum 
requirement. More specihcally, aeA ^iSr i — Vk where 

rjk denotes the target share of the airtime for ISP k. 

Accordingly, to formulate such optimization problem, the 
feasibility region of t“ (or xf) is required. Thus, we study 
the behavior of a single STA -which is using EDCA- with a 
three-dimensional Markov chain. With the aid of the proposed 
Markov model, we can learn how to implement or control 
t“ in terms of EDCA parameters. As a result, based on the 
established relationship between r“ and EDCA parameters, 
we would be able to analyze its feasibility region and also 
design an algorithm to control EDCA parameters to approach 
the optimal r“. 


III. Optimization Problem 

In this section, we present the STA-AP association and 
airtime control optimization problem based on the system 
model introduced in Section |II] We first use a Markov chain 
model to model the EDCA protocol. With this model, we 
can analyze the feasibility region of rf for each STA. Then, 
the optimization problem is formed and solved by applying 
monomial approximation and geometric programming. 
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Fig. 4: Three dimensional Markov chain model for one STA 

A. Markov Chain Model for EDCA 

To study the feasibility region of rf and its implementation 
in the EDCA MAC protocol, we establish a three-dimensional 
Markov chain. Such model enables us to estimate the transmis¬ 
sion probability of the STAs in a WLAN employing the EDCA 
protocol. Different from the previous works, which discuss a 
single parameter in the EDCA protocol (e.g., fTS)), we aim at 
developing a general mathematical model that can show the 
influences of all the EDCA parameters on the transmission 
probability. 

Eigure shows the three-dimensional Markov chain model 
to describe the behavior of a single STA, which uses the EDCA 
MAC protocol. Since STA transmissions are only coupled 
within a BSS, we consider a single BSS where an AP and its 
associated STAs reside in. Thus, in this subsection without loss 
of generality, we remove the index of a (e.g., rf is replaced 
by Ti) to keep the notation simple. 

Our proposed three-dimensional Markov chain model is an 


extension to the presented model in pT) . To be able to further 
control the competition among STAs and thus the collision 
probability in a BSS, we introduce two new MAC parameters 
(i.e., Qi and Li for STA i) for each STA. Similar to p4) , after 
a successful transmission or a packet drop, the STA has to 
flip a biased coin with successful probability qi to enter the 
backoff process. Whenever the STA fails to enter the backoff 
process, it needs to wait for a period of long inter-frame space 
(i.e., Li time-slots) before another try. These two variables are 
incorporated in the proposed Markov chain. 

It should be noted that time in this Markov chain is slotted 
and the interval between any two adjacent states is a general 
time-slot. Here, a general time-slot can be an idle time-slot 6, 
a successful transmission duration (Tg) or a collision duration 
(Tc). In this Markov chain, the triple {s(f), b{t), u(f)} denotes 
the state of the STA at time t. More specifically, s{t) represents 
the backoff/retransmission stage of STA, where —2 < s{t) < 
rrii + hi (s{t) = —2, —1 will be explained later). Furthermore, 
b{t) denotes the backoff counter, which takes a value between 
0 and the contention window of the current backoff stage. 
Furthermore, v{t) represents the remaining time of a frozen 
period or the waiting time to the next trial when the STA 
already failed to enter the backoff process. 

The detailed transition probabilities between states of the 
Markov chain and the stationary condition of the Markov 
chain can be found in the Appendix. With such stationary 
conditions, the transmission probability t of each STA is 
derived as function of the detailed EDCA parameters as in 
0, where pi is the collision probability faced by STA i as 
defined in (|26ll. 


B. Feasibility Region Analysis 

In this subsection, we study the feasibility region of rf 
in order to complete the optimization problem formulation. 


Based on (13i, the adjustable variables to tune t“ are to“. 


hf, Lf, qf, Af, and Wf, where 0 < ?' < mf + hf. Recall 


that 0 < < 1, 0 < L“, and 0 < Wfi. Furthermore, since 


A, = AIFS[STA i] - 1, and AIFS > SIFS in EDCA , we have 
Ai > SIFS. Assuming that SIFS is equal to one time-slot, we 
have A“ > 1 without loss of generality. 

J) Lower Bound of rf: Let r“ be the lower bound of rf. 

0 when —>■ oo, or qf —> 0, 


From (13 1 , it is clear that rf 
or A“ —> c». Consequently, 


r“ = 0 


(14) 


2) Upper Bound of rf: Let rf be the upper bound of rf. 
To calculate r“, we need to make the denominator of r“ in 
































0 as small as possible. Thus, we set Lf = 0 (or qf = 1) 
and j = 0, which make the first and fourth terms in the 
denominator zero. Then, by setting Ai to its lower bound (i.e., 
= 1), we can minimize the second term in the denominator 
of t“. Consequently, at Lf = 0, ^ = 0, and Af = 1, 


1 + 
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(15) 


Thus, by applying successive transformation strategies, we 
will try to convert the original problem into a series of 
standard GP problems that can be solved to reach an optimal 
solution. First, we introduce three auxiliary variables, = 
nj'g 5 (l + xf,) — t', \/a G A, uf = 1 — pf, Vi e S,Va G A, 
and tf = 1 + xf, Vi G S,Va G A. Then, the optimization 
problem in ([T7| can be transformed into 


From (15 1 , it is clear that the upper bound of t“ can be 
achieved when mf + hf oo. Thus, ff can be written as 
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C. Optimization Problem 

Taking into account the feasibility region of r“ derived 
in Section III-B[ here, we can mathematically present the 
transmission probability optimization problem. The objective 
is to maximize the overall network throughput, while dis¬ 
tributing access airtime among different ISPs according to 
their reservations. More specifically, the optimization can be 
formulated as 

E X^T^t 
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(18f) 

(18g) 


where T = [tf] (tf > 1), U = [uf] (0 < uf < 1), Y = [y“] 
(?/“ > 0). Nevertheless, the transformed problem is not still 




> rjk, Vk G JC (17b) in a GP form. One reason is that the objective function (18a I 
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where X = [xf] (xf > 0) and P = [pf] (0 < pf < 1). Let 
us recall that xf = jzfpc and pf = 1 — (1 ~ xf). It 

should be noted that the optimization problem is alternatively 
formulated with respect to xf instead of t“, since it will prove 
useful to solve the problem. 

In the preceding optimization problem, the objective func¬ 
tion in (|17a|l represents the overall network throughput (i.e.. 


is not a posynomial because of negative multiplicative coeffi¬ 
cients. To deal with such problem, first, we equivalently substi¬ 
tute the objective function by a^A ~xfxft{y°‘)~^ + M 

where M is a sufficiently large positive constant. Adding 
M makes sure that the objective function is always positive. 
Then, we introduce an additional auxiliary variable xq > 0. 
By minimizing xq and guaranteeing constraint Cll in (|I9 1 , 
we can effectively minimize the objective function in ( |18a i. 
Consequently, 


min xq, subject to, 
X,T,U.P,Y,xo 


(19a) 




Tf) based on (10 1 . Furthermore, constraints in 


( |17b[ ) guarantee the minimum airtime reservations for all ISPs 
(i-e-. > Vk) based on ([T^. This set of 

constraints enable controlling ISPs’ share of access airtime 
regardless of their number of STAs. To guarantee that the 
provided solution fits into the feasibility region, in constraint 
(17c I, t“ is limited to its upper bound. Then, the equality 
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constraint (17di establishes the relationship between the upper 
bound ff and pf based on ([T6|). Finally, constraint (17e i estab¬ 


lishes the relationship between pf and all xf, Vi' G S, i' i 
according to ( p^ . 

The formulated problem is non-convex and thus intractable 
to solve. However, it potentially looks like an extension of 


C16 : 


i' 

=l,ViGS,aGA 


Geometric Programming (GP) (defined in Section VII-B i. 


I + X 

In the preceding optimization problem, pf is replaced by l—uf 
based on constraint in (|18g|l. The optimization problem in 
0 belongs to the class of complementary GP problems 













































Algorithm 1 : GP-based Association Control Algorithm 

Initialize tf, pf, uf, for all i S 5, a S A, ya for all 

A, xq\ 

Record the current system state as Z = (X, T, U, P, Y, ccq); 

repeat 

Compute the ratio a of each monomial term in the 
denominator of Cll, C12, C13, C14, and C16 according 
to ( |30| ), at the current system state Z; 

Apply monomial approximation to the denominators 
mentioned above according to ( |29] l; 

Solve the resulting GP problem using cvx; 

Update the current system state Z = {X,T,U, P,Y, xq); 
until all xf converge. 

Compute the optimal transmission probabilities r*“ = 


that allow upper bound constraints on the ratio between two 
posynomials and equality constraints on the ratio between a 
monomial and a posynomial p5) , By approximating 

the posynomials in the denominator of such constraints, a 
complementary GP can be turned into a standard form of 
GP. Consequently, the optimal solution can be achieved by 
iteratively applying monomial approximations and solving 
a series of GPs. The arithmetic-geometric mean inequality 
can be used to approximate a posynomial with a monomial. 
The details of such monomial approximation are provided in 
Section IVII-BI 

Accordingly, we propose an iterative algorithm to reach 
to an optimal solution of the transmission probability opti¬ 
mization problem. In each iteration, monomial approximations 
are applied to the denominator of Cll, Cl2, C13, Cl4, and 
C16. Then, the resulting GP can be solved for instance by 
using a standard interior-point algorithm. More specifically, 
Algorithm[T]presents different steps need to be performed until 
convergence. 


D. Asymptotic Complexity Analysis 

In this subsection, we study the asymptotic complexity and 
scalability of Algorithm in terms of the number of STAs 
(i.e.. Ns) and the number of APs (i.e.. No). 

In Algorithm for each iteration, the computational com¬ 
plexity is incurred by applying monomial approximations 
and solving the resulting GP problem. Suppose Cma and 
Cgp denote the required computing efforts for monomial 
approximations and solving GP problem in each iteration. 

More specifically, first in each iteration, the denominators 
of constraints Cll, C12, C13, C14, and C16 need to be 


approximated according to (29 1 . The required computational 
complexity for these monomial approximations is proportional 
to the number of monomial terms in the denominators of 
constraints Cll, C12, C13, C14 and C16. Consequently, the 
complexity of monomial approximations totals to 


Cma =OiNaNs) + C(iV„iV«) + C(iV,iVf) 

+ 0{NaNs) + 0{NM = 0{NaN^) 


( 20 ) 


Subsequently, the resulting GP problem needs to be solved 


by transforming it to a convex problem. It is reported in p7| 
that the worst-case computational complexity of this approach 
is 0{pn?), where n is the number of variables, and p is the 
total number of terms in all monomials and posynomials in the 
objective function and constraints. The number of variables of 
the optimization problem in ( [T^ is 

n = 47V„iV, + 27V„ + 1 = C(X,7V,) (21) 


Furthermore, the total number of terms in all monomials and 
posynomials can be counted as 


P —1 + {4:NaNs -|- 1) -f Na{Ns + 2) 

+ + NkNaiNs + 1)) 

k 

+ 7NaNs + NaN^s + 

=0{NaN^s) ( 22 ) 


where each term in ( [22] l is respectively the number of terms in 
all the monomials and posynomials in the objective function 
and constraints C11-C16. Based on ( |2T] l and ( |22| ), to solve the 
GP problem, the complexity is 

Cgp = 0{pn^) = 0{NtNl) (23) 


Thus, based on ( |20[ ) and ( |23| l, per-iteration asymptotic com¬ 
plexity of Algorithm becomes 

Cai^ iteration — Cma + Cgp = OiN^Nl) (24) 


We now turn to studying the number of iterations (denoted 
by A) required for Algorithm to converge. Figures 5a and 
l5bl illustrate numerical results on A versus number of APs 
(i.e.. No) and average number of STAs per AP (i.e., Amean), 
respectively. The simulation setup used in these figures is the 
same as the one presented in Section |V] In Figure 5a 




is set to so that the expected total number of STAs in 
the WLAN stays fixed. Figure shows that A increases with 
Na when Na is small, and then, it becomes steady for a larger 
range of Na- In Figure]^ Na is fixed equal to 4 and Amean is 
varying. Figure shows that A grows linearly over a typical 
range of Amean- It should be noted that the fluctuations in A 
are caused by the randomness in the number of STAs per AP, 
which follows a Poisson distribution with mean of Amean- 
Therefore, based on ( |24| ) and the numerical results in Fig¬ 
ures and it is shown that the overall complexity of the 
proposed algorithm only grows polynomially with the number 
of STAs and APs, compared to an exponential complexity 
required by direct search methods. 


IV. Implementation Details 
The optimal transmission probability of STA i at AP a (i.e., 
T*“) can be obtained from Algorithm But, in the EDCA 
protocol, the transmission probabilities of STAs are not di¬ 
rectly controllable. Instead, what we can control are the MAC- 
layer parameters (e.g., IFmin and AIFS) to achieve the optimal 
performance. In this section, we first verify the accuracy of the 
relationship between r“ and the EDCA parameters provided 
in ([T3 ]i and also the validity of the Markov chain model. Then, 













(a) Number of iterations vs. number of APs (Na) Fig- 6: Accuracy of controlling different parameters 



(b) Number of iterations vs. STA density (Amean) 

Fig. 5: Number of iterations required for Algorithmto converge 


we develop an algorithm to adjust EDCA parameters aiming 
to implement in the EDCA protocol. 

We first investigate the achievable accuracy by controlling 
different parameters through an example. In this example, we 
consider one BSS with 6 STAs. The transmission probabilities 
of STA 2 to STA 6 are fixed (e.g., T 2 = T 3 = T 4 = = tq = 

0.005), while the transmission probability of STA 1 (i.e., ri) 
is varied from 0.005 to 0.1. 

Eigure compares the throughput of STA 1 analytically 
derived from ( [T0| i and numerically measured using an EDCA 
simulator. In the numerical results, different EDCA parameters 
are separately adjusted to tune ti to the desired value. Erom 
analytical result in Eigure it is clear that the throughput 
of STA 1 increases with its transmission probability ti. 
Regarding the numerical results, the following conclusions can 
be drawn. 

Erom Eigure FEmin.i is the best single parameter to be 
controlled aiming to adjust to a desired value. Ai and Li 
provide a good match and thus can be good choices to tune 
Ti. Eurthermore, rrii and hi control is not adequately effective, 
especially for hi. Also, as we can observe, the achievable 
ranges of r by controlling each parameter (i.e., Ai, W^nin.i and 
Li) are limited. More specifically, when Ai is being adjusted. 


Algorithm 2 : MAC Parameter Control Algorithm 


For each STA-AP pair {i, a): 
Initialization: 


= 15, A“ = 6, qf = 0.5, Lf = 100, mf = 6, 


Set 

= 6 ; ^ 

Compute j from (131, while other variables are fixed; 
Round i to an integer; 
if i < 6 then 

Set lT“i„ j = 0; Compute 


from (13 1 , while other 


variables are fixed; 

if < 0 then 

Set Lf = 0; Compute Af from while other 
variables are fixed; Round Af to an integer; 

if A“ < 1 then _ 

Set A“ = 1; Compute mf 


from (13 I, while other 


variables are fixed; Round mf to an integer; 

if mf < 0 then 


Set mf = 0; Compute hf from (131, while other 
variables are fixed; Round hf to an integer; 

li hf <0 then 

Set /i“=0; 

Start transmission using EDCA protocol with parameters 
Af, qf, Lf, mf, hf}. 


we can learn from Eigure that the simulation result matches 
with the analytical result when ri is in the range [0.005,0.15). 
When CWmin is being controlled, the situation is similar, 
but for a broader range of ([0.005, 0.065)). Eigure also 
suggests that controlling L is effective when ti is in the range 
[0.005, 0.035). Thus, in summary, to achieve a boarder feasible 
range of r, we need to jointly adjust these parameters. 

Here, we aim to develop a MAC parameter control algo¬ 
rithm to achieve the optimal transmission probability obtained 
from the optimization problem in Section |I^ In Section 
the feasibility region of transmission probability is calculated 
assuming that all the 5 MAC parameters (i.e., Wmin, A, q or 
L, m, h) can be freely varied in their feasible ranges. Thus, a 
control algorithm is required which can be able to adjust all the 
5 parameters simultaneously. Thus, an algorithm is proposed 
taking into account the achievable accuracy by controlling 
different parameters. Algorithm shows the complete control 
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Modulation 

PEC Rate 

Data Rate (Mbps) 

SNR (dB) 

BPSK 

1/2 

6 

[5,8) 

BPSK 

3/4 

9 

[8,10) 

QPSK 

1/2 

12 

[10,13) 

QPSK 

3/4 

18 

[13,16) 

16QAM 

1/2 

24 

[16,19) 

16QAM 

3/4 

36 

[19,22) 

64QAM 

2/3 

48 

[22,25) 

64QAM 

3/4 

54 

[25,00) 


TABLE I: IEEE 802.11a adaptive modulation and coding scheme 
and the SNR ranges used in the numerical results 


procedure. 

The achievable throughput for STA 1 in Algorithmic is also 
plotted in Figure As expected, we can learn that Algorithm 
1C can approach the desired throughput closely for the total 
range of ti (i.e., the upper bound ti computed from ([T6|)). 
Thus, Algorithm|Ccan improve the control accuracy compared 
with controlling the parameters separately. 


time-slot 5 

9 gs 

Propagation Delay 7 

1 gs 

Ttxop 

1 ms 

SIPS 

10 gs 

ACK 

40 gs 


TABLE II: IEEE 802.1 le MAC parameters used in the numerical 

results 
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V. Numerical Results 


In this section, we present numerical results to evaluate 
the performance of the proposed STA-AP association and 
airtime control algorithm, and also the MAC parameter control 
algorithm. More specifically, the performance of our GP-based 
association scheme is compared with the Max-SNR scheme in 
terms of throughput and fairness. We implemented a simulator 
for EDCA (including access probability q and inter-frame 
space L described in Section VII-A| i in Matlab to measure 
the achieved throughput numerically. Furthermore, we used 
cvx to solve the GP problems in the association algorithm. 

We consider a network in which 4 APs are deployed in a 
10 X 10 area. More specifically, the APs are placed at the 
centers of four different 5x5 m? grids to provide seamless 
coverage. To eliminate interference between the transmission 
of different APs, four non-overlapping 20 MHz channels are 
assigned to four APs. The STAs are distributed in the entire 
area according to the two-dimensional Poisson point process 
(PPP). 

The wireless channel model includes path loss and small- 
scale fading. Generally, the channel gain can be expressed as 
h = where d is the distance between a STA and an 

AP, a > 2 is the path loss exponent, A is a constant dependent 
on the frequency and transmitter/receiver antenna gain, and h' 
represents the small scale fading component. In the numerical 
results, we set a = 3 and A = \. Furthermore, h! is randomly 
generated according to the Rayleigh distribution assuming 
= 1- The received SNR at STA i is equal to 
where P is the transmission power, gf = \hf\'^ is the channel 
power gain from STA i to AP a, and cr^ is the power of noise. 
In all the numerical results, Pja^ is assumed fixed and set to 
lOdB. 


To determine the transmission rate of each STA-AP pair, 
the 802.1 la physical layer model is used. More specifically, to 
guarantee a maximum packet error rate, adaptive modulation 


Fig. 7: Throughput vs. STA density, homogeneous distribution 


Throughput vs Density: GP (Nonhomogeneous Distribution) 



Throughput vs Density: Max SNR (Nonhomogeneous Distribution) 



Fig. 8: Throughput vs. STA density, non-homogeneous distribution 

and coding is used based on the received SNR. Table |I] shows 
the achievable transmission rates and adaptive modulation and 
coding schemes standardized in IEEE 802.11a and the SNR 
range for each scheme used in the simulations. 

The MAC layer parameters used in our simulations are 
summarized in Table Moreover, the target airtime share 
for each ISP k (i.e., rjk) is set equal to the number of APs 
divided by the number of ISPs. In other words, we assume 
that the ISPs have the same minimum airtime reservation and 
share the total airtime in a fair manner. In the Markov chain, 
the average frozen time N is approximated by Ttxop/^- 

A. Ejfects of STA Distribution 

Here, we investigate the impact of STA distribution on 
the fairness and throughput achieved by the two association 
algorithms. More specifically, we set up two examples consid¬ 
ering homogeneous and non-homogeneous STA distributions. 
In both examples, the STAs are randomly associated to the 
two ISPs with equal probabilities. 


























































ISP, load (p,) 

s Distribution Balance: Max SNR 



ISP, load (p,) 


Fig. 9: Throughput vs. ISP 1 load 

1) Homogeneous Distribution: In this example, the STAs 
are distributed in the square space according to a homogeneous 
two-dimensional PPR Accordingly, the number of STAs in 
each grid (where an AP centered at) follows a Poisson distri¬ 
bution with mean Amean, which represents the average number 
of STAs per AP. Figure |7] shows that GP-based association 
scheme improves both fairness and total throughput compared 
with Max-SNR scheme. 

2 ) Non-homogeneous Distribution: Let consider that STAs 
are distributed according to a non-homogeneous PPP More 
specifically, the number of STAs in a grid (where AP a 
centered at) follows a Poisson distribution with mean A^. 
Here, Aa is randomly generated between 0 and Amean- With 
such STA distribution. Figure]^ shows the achieved throughput 
of two ISPs versus Amean- As expected, the performance 
gap is significantly larger with non-homogeneous distribution 
compared to the homogeneous case. Max-SNR can hardly 
guarantee the fairness between different ISPs. While for the 
same distribution, GP-based STA-AP association can man¬ 
age to keep the balance between two ISPs. Furthermore, it 
can offer improvement over Max-SNR in terms of the total 
throughput due to the load-balancing among APs. 


B. Effects of STA Density and ISP Load 

Let define pi (also referred to as ISP 1 load) as the ratio of 
number of STAs serving by ISP 1 to the total number of STAs 
in the network. Here, the performance of the two association 
approaches are compared under different STA density and ISP 
load. 

Assuming a homogeneous STA distribution with Amean = 3, 
Figure demonstrates the achieved throughput of two ISPs 
versus different values of pi for both GP-based and Max- 
SNR association schemes. By Max-SNR association, it is 
shown that throughput of ISP 1 grows linearly with pi, while 
the achieved throughput of ISP 2 is decreasing. But, GP- 
based association can fairly distribute the airtime between ISPs 
regardless of their ISP loads and thus keep a perfect balance 
between the achieved throughput of the two ISPs. 

Figure 10 shows the total throughput achieved by the two 
association algorithms versus Amean for a homogeneous STA 
distribution. For a fixed pi, the total throughput by both 



Fig. 10: Total throughput vs. STA density for different pi 



Fig. 11: Fairness vs. STA density for different pi 


algorithms increases with the STA density (i.e., Amean)- But, 
the throughput increase rate is decreasing with Amean- This 
is because the wireless channel is underutilized when the 
STA density is low. Thus, the increase in the STA density 
will improve the total throughput. But, when the STA density 
is large, increasing the STA density further will result in a 
higher collision probability, and hence, slow down the total 
throughput improvement. For any fixed pi, it is shown that GP- 
based association significantly improves the total throughput 
compared with the Max-SNR association. 


Figure 11 measures the fairness by employing the Jain’s 
fairness index in 


(Sfcg/C 

1^1 ^keK 


(25) 


where agA '^i achieved throughput for all 

the STAs of ISP k. From Figure[^ it is clear that the proposed 
GP-based association approach can always guarantee perfect 
fairness between the ISPs regardless of the STA density or pi. 
The achieved fairness level by Max-SNR association is always 
worse than GP-based, especially when the STA load is highly 
unbalanced between ISPs (i.e., pi is not close to 0.5). 








































( 26 ) 


VI. Conclusion 

This paper considers the STA-AP association and airtime 
control in virtualized 802.11 networks aiming to provide 
fairness guarantees among ISPs despite the number of STAs 
per ISP. First, a three-dimensional Markov chain is developed 
to model a generalized 802.1 le EDCA protocol. This model 
establishes the relationship between the transmission proba¬ 
bility of each STA and the detailed parameters in the MAC 
protocol. Based on this relationship, the feasibility region of 
the transmission probabilities are derived. Subsequently, an 
optimization problem is formulated which can maximize the 
network throughput, while guaranteeing the fairness between 
different ISPs. The implementation of the optimal transmission 
probabilities obtained by successive geometric programming 
is discussed by controlling the MAC parameters. Extensive 
numerical results confirm that the Markov chain can accurately 
describe the MAC protocol. Eurthermore, it is verified that the 
proposed association algorithm can improve the throughput 
and provide fairness guarantees in virtualized 802.11 WLANs. 

VII. Appendix 
A. Markov Chain Model 


probabilities of all STAs (i.e., ) in the BSS, we have 

P* = 1 - n (1 - A/). 

i' 

Eor the backoff process, we use Eigure [T^ to explain the pos¬ 
sible transitions at the backoff/retransmission stage j. When 
STA i reaches the state {j — 1,0,0}, it will start a transmission. 
Then, if other STAs happen to access the channel at the same 
time, STA i will experience a collision. The probability that 
STA i encounters a collision is equal to pi. In a case of 
collision, the STA will enter the next backoff/retransmission 
stage (i.e., j), pick a random backoff counter b between 0 and 
Wj^i, and go to the state {j, b, 0}. When the STA reaches the 
state {j, b, 0 }, it will count down to {j,b—l, 0 } if the channel 
is idle. Otherwise, if the channel is sensed busy (i.e., other 
STAs are transmitting in the channel), STA i will freeze its 
backoff counter and move to the state {j, b,N + Ai}. In the 
frozen period, the STA will count down v{t) by 1 each time- 
slot until it gets to the state {j, b,Ai}. Then, if the channel 
is still idle (with probability (1 — pi)), the STA will continue 
counting down to {j, b, 1 } and then go to state {j, 6 — 1 , 0 }. 
Otherwise, if the channel becomes busy when STA i is in the 
state {j, 6, d} where 1 < d < A^, it has to freeze its backoff 
counter again and go back to the state {j, b,N + A^}. 


Here, we detail the Markov chain model for EDCA intro¬ 
duced in Section III-A Eirstly, the transition rules among 
the states of the Markov chain are given according to the 
EDCA protocol. Then, the transition probability between any 
two adjacent states are summarized. Einally, by solving the 
stationary point of the Markov chain, the relationship between 
transmission probability and the detailed EDCA parameters 
are derived. 


1) Transition rules: Eirst, we describe the transitions 
among the states in the proposed Markov chain according to 
the EDCA protocol. As shown in Eigure after a successful 
transmission or a packet drop, the STA will flip a biased 
coin with successful probability qi. If the coin appears to be 
tail (with probability 1 — qi), the STA will go to the state 
{—2,0, Li — 1}, and try again after Li time-slots. Otherwise, 
if the coin is head (with probability qi), the STA will move 
to the state { —l,0,Ai}, where Ai = AIES[STA i] — 1. If 
the channel stays idle in the following Ai + 1 time-slots, the 
STA will move on to the backoff/retransmission stage 0 and 
uniformly pick a random backoff time from [0, Wo,i]. 

Otherwise, if the STA senses the channel busy in a state 
{ — 1 , 0 , d} where 0 < d < Ai, it will freeze and move to 
the state { —1,0, A^ + Ai}, where N is the average frozen 
time (The frozen time is equal to either Tg or and thus N 
can be approximated to T as discussed before). Afterwards, 
it will wait until channel becomes idle again, by counting 
down v{t) from N + Ai to A,;. When the STA gets to 
the state { —l,0,Ai}, it has to keep repeating the explained 
process until it enters the backoff process. We assume that 
the probability that the channel is observed busy by STA i is 
a constant pi in a general time-slot. Given the transmission 


Einally, for states {j, 0,0}, 0 < } < rui + hi, the STA 
will initiate a transmission. If no other STA transmits at the 
same time (with probability (1 — pi)), the STA will have a 
successful transmission of duration Tg. Otherwise, if any other 
STA simultaneously starts a transmission (with probability pi), 
STA i will experience a collision of duration T),. After the 
collision, if 0 < j < rrii + hi — 1, the STA will go to 
the next backoff/retransmission stage. But, if j = rrii + hi 
(which means that the maximum backoff/retransmission limit 
is reached), the STA has to drop this packet. Then, for a new 
packet transmission, it has to flip a biased coin and start over 
from the state {— 2 , 0 , Li — 1 } or {— 1 , 0 , A^j according to the 
result of the coin toss. The mathematical representations of 
the transition probabilities are summarized in Section [VII-A2| 

2) Transition Probabilities of the Markov Chain: Given the 
transition rule between the states described in above, we have 
the following transition probabilities. 

1) Eor states {—2,0,d}, 0 < d < — 1, i.e., after getting a 

tail in a coin toss, 

rP{(- 2 , 0 , d - l)|(- 2 , 0 , d)} = 1 , 1 < d < L, - 1 
<^P{(-2,0,L,-1)|(-2,0,0)} = 1-(Z, 
[p{(-l,0,A,)|(-2,0,0)} = g, 

2) Eor states {—1,0, d}, 0 < d < N + Ai, i.e., before entering 
the backoff process, 

{ P{(-l,0,d- l)|(-l,0,d)} = 1, A,-f 1 < d< A^-f A, 
P{(-l,0,d- l)|(-l,0,d)} = l-p^, 1 < d< A, 
P{(-1,0,W +A,)|(-l,0,d)} =K, 0 < d< A, 
P{(0,6,0)|(-l,0,0)} = pj^, 0<6<lVo,. 







3) For states {j, b,d}, 0 < j < rrii + hi, i.e., in the backoff 
process, 

a) If 6 = 0,c? = 0, i.e., before a successful transmission or 
a collision, 

'P{(-l,0,A,)|(j, 0,0)} = (l-p,)ft, 

0 < J < rrii -\- hi — \ 

P{(- 2 , 0 , U - l)|(j, 0 , 0 )} = (1 - pi){l - Qi), 

0 < j < rui + hi — 1 

' p{(j + i,&,0)|O',0,0)}=p,pp^^, 

0 < j < rrii + hi — 1, 0 <b < VFj+i,i 
P{{-2,0,Li - l)\{mi + /ii,0,0)} = l-qi 
P{{-l,0,A^)\{mi + = qi 

b) If d = 0,1 < 6 < Wj^i, i.e, in the count down process, 

f P{{j, b - 1, 0)|(j, 6, 0)} = 1 - p„ 1 < & < Wj,, 
\p{{j,b,N + A,)\{j,b,0)} =pi, l<b< Wj^i 

4) If I < d < N + Ai, I < b < Wj^i, i.e., during a frozen 
period or the AIFS after it, 

f P{{j,b,d- l)\{j,b,d)} = 1, A, + l<d<N + Ai 
I P{{j,b,N + A,)\{j,b,d)} =pi, l<d<Ai 
I P{U:b,d- l)\{j,b,d)} = l-pi, 2<d<A, 

[p{U:b- 1,0)10', 6,1)} = 1 - Pi 


By exploiting the transition structure within a frozen process 
and the post AIFS (i.e., the states {j, b, d}, 0 < d < N + Ai), 
we have. 


bj,b,d — 
bj,h,d — 


(1 0 < j < + /li - 1,1 < d < Ai - 1 

(1 -p-)Ai f)<j<m^ + hi,A^<d<N + Ai 


For states {—2, 0, d}, by looking at the incoming and outgoing 
probability, we can have 


b-2,0,d 


- ~bo 0,0, 0 < d < Li — 1 

Qi 


For states {—1,0, d}, by analyzing the structure of transition 
between the states {—1, 0, d}, 0 < d < N + Ai and the state 
{0,0, 0}, we can have 


^- 1 . 0 ,d 
^- 1 , 0 ,d 


7 j-0 < d < Ai 

_ (\ _n 

— 7z -F4~+i—^ 0,0,0) Ai + l<d<A^ + Ai 


To this point, all the stationary probabilities in the Markov 
chain are represented in terms of bo^o.o- Then, &o,o,o can be de¬ 
rived from the normalization condition in p7| ). Consequently, 
having all stationary probabilities &j,b,d, the transmission prob¬ 
ability of STA i can be calculated as in To derive 0, 
the contention window Wj^i is set equal to 


3) Stationary Condition of the Markov Chain: Based on the 
transition rule above and the structure of the Markov chain, 
as in nn and 1 ^ , we have, 

bi,o,o = Plbo,o,o, 0< j < rrii + hi 
Wa + l-b 

bj b,o “ jxr I t 'bj op, 0 ^ J ^ rrii hi, f A b f. Wj i 

Wj,i + 1 


pp _ ibkmin.iS-t, 0<j<mi 

\ , mi + l<j <mi +hi 

according to the exponential backoff rules in 802.1 le. From 
the definition of transmission probability, can be calculated 
by summing all the stationary probabilities of the states in 
which the STA will initiate a transmission (i.e., bjpp, Vj € 


















Li —1 N-\-Ai rrii+Li mi-\-hiWj^i N-\-Ai 

1 = ^-2,o.d + b-lfi,d + bo,ofi+ Y Y Y 

d—0 d—O j—O j—O b—l d—O 



1 +PiN 1 - (1 

Pz (l-p*)^>+i 


1 -p: 


mi-\-hi-\-l 


1 -p* 


1 + Np^ 

2(1 


mi-\-hi 

Y 


3=0 


bo,0,0 


(27) 


{0,1,..., TOi + hi}). The expression in (13i establishes the 
relationship between and all the detailed parameters for STA 
i in the EDCA protocol, which will be useful to characterize 
the feasibility region of rf in the optimization and design an 
algorithm to control the transmission probabilities of the STAs. 


B. Geometric Programming and Monomial Approximation 

An optimization problem is called geometric programming 
if it follows the following form, 

min/o(x), subject to 

X 

fi{x) <1, i = 1,... ,ni 

gi{x) = 1, J = 1,... ,712 

where fo, ■ ■ ■, fm posynomials and gi,, gn^ are mono¬ 
mials. In the context of geometric programming, a monomial 
function / of x = (xi, X 2 ,..., x„) is defined as, 

fix) = CX“1X2" . . 


where c > 0 and Oi € TZ. Furthermore, a posynomial is defined 
as the summation of multiple monomials, i.e., 

5(^) = 

The basic idea of monomial approximation is as follows: 
consider a posynomial function g{x) = fk{^) with /fc(x) 

being the monomial terms. By the arithmetic-geometric mean 
inequality, we have 



Otki^o) 


(29) 


where the parameters a/j(xo) can be obtained by computing 


ctkixo) 


fkjxo) 

gi^o) ’ 


Vfc 


(30) 


where Xq > 0 is a fixed point (e.g., the optimal solution from 
the last round of optimization). It is proved that g{x) is the 
best local monomial approximation of p(x) near Xq 
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